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Abstract. A scenario of heavy resonances, called massive Hagedorn states, is 
proposed which exhibits a fast [t rj 1 fm/c) chemical equilibration of (strange) baryons 
and anti-baryons at the QCD critical temperature T c . For relativistic heavy ion 
collisions this scenario predicts that hadronization is followed by a brief expansion 
phase during which the equilibration rate is higher than the expansion rate, so that 
baryons and antibaryons reach chemical equilibrium before chemical freeze-out occurs. 
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1. Introduction 



The enhancement of (multi-)strange (anti-)baryons has been predicted as a potential 
proof for the existence of the quark gluon plasma (QGP) in relativistic heavy ion 
collisions pp. Hadro-chemically saturated multiplicities of (anti-)hyperons have been 
experimentally discovered in central collisions in Pb+Pb experiments at CERN-SPS 
and in Au+Au experiments at Brookhaven, see Ref. [2|. 

At SPS energies, the strong increase in the yield of antiprotons |3] and antihyperons 
|H Ej has been explained by a "clustering" of mesons of the following type: 

riiir + n 2 K <-> Y+p. (1) 

This channel is very efficient to chemically equilibrate (strange) antibaryons in a baryon- 
dense fireball as is expected at the SPS energies. The equilibration time scale is 
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inversely proportional to the density of baryons and to the annihilation cross-section, 

1. e., (Ty)^ 1 ) ~ \/{a B YV B Yn B ). Assuming that various cross-sections are related as 
given roughly by the additive quark model, &BY->nir+n Y K ~ o~Np^nn, (for a theoretical 
model, see Ref. |6j) and that the density is about 1 to 2 times larger than the normal 
nuclear matter density po, one finds equilibration times for antibaryons t eq ~ 1-3 fm/c 
0j. This leads to microscopically calculated antihyperon yields which are consistent 
with the fitted chemical freeze-out parameters 0. 

The multi-meson reactions at RHIC cannot account for an efficient, direct baryon- 
antibaryon pair production within the standard hadron resonance gas model. The 
observed "chemically saturated" (strange) (anti-)baryon yields |H1 IH Ej cannot be 
obtained from hadron cascades. Antibaryons are much more abundant at RHIC than 
at SPS, and equilibration times are ~ 10 fm/c at a temperature of about 170 MeV 
IHEIIZI. Nevertheless, the RHIC hadron multiplicities suggest that chemical equilibrium 
is reached at about 170 MeV ~ T c . Thus, it has been suggested that various hadron 
species are "born in equilibrium" 

Here, a scenario of dynamical equilibration is developed, overcoming some 
shortcomings (see Sec. EJ) of Ref. [H] (hereafter referred to as BSW). According to 
BSW, the reactions in Eq. (JTJ alone lead to chemical equilibration at RHIC. Sec. |3] 
presents a sufficiently fast (~ 1 fm/c) chemical equilibration mechanism for both of 
(strange) baryons and anti-baryons. It relies on abundant (at T c ) massive resonances, 
called Hagedorn states (HS), which are short lived. (Strange) baryon and antibaryon 
production proceeds along the reaction 

(tiitt + n 2 K + n 3 K ^) HS ^B + B + X. (2) 

Here X represents all possible multi-hadron states. Such reactions are an important 
generalization of Eq. (£Q). The equilibration time can be estimated by the branching 
ratios of HS decays into B + B + X using a microcanonical statistical model. We then 
address the production of strange (anti-)baryons, especially the rare Q (sss) state. The 
role of HS for chemical equilibration at nonzero net baryon density (e.g., at SPS and 
AGS energies) is briefly discussed. 

2. Overpopulated hadron densities 

Typical equilibrium baryon/antibaryon densities in a net baryon-free hadronic system 
at T pa 170 — 180 MeV are n B = n'g ~ 0.04 fm~ 3 (neglecting eigenvolume effects). 
With annihilation cross sections (o~v) ~ 30 mb, the equilibration time due to reactions 
as in Eq. (jj) is t eq « 10 fm/c. This cannot explain the apparent chemical equilibration 
of baryon and antibaryons jUEUIJ. It was suggested in Ref. [§] that such multi- meson 
collisions can still lead to quick chemical equilibration, in close vicinity to the phase 
transition: by comparing a hadron gas model equation of state with the lattice results, 
it was speculated that there exists a state of extra-large particle density around T c , which 
is effectively overpopulated with pions and kaons. This overpopulation then could drive 
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the baryon and antibaryon pair production on a very short time scale. In essence, the 
total production of a specific type of baryon (e.g., the Q as the most exotic one) is 
attributed |Hj to the gain term in the following master equation [U El El Ej : 

^-n n = -"loss" + "gain" = -(a nS v nS ) x (n s n n - ]T M {num) {n^{n K ) n A . (3) 

Here the mass-action factor reads M( nij „ 2 ) = ti^tiq /(n^) ni (n^) n2 . Notice that, 
in accord with Eq. Q, the density of Qs tends to approach during chemical 

equilibration. If the equilibrium is maintained, i.e., when the "loss" and "gain" reactions 
have equal rates, saturation does not change. In this case, the production (as well as 
the annihilation) rate are given by ~ {crnB v ciB) n Bi an d the corresponding chemical 
equilibration time is Tq ~ 1/Tq. This gives > 10 fm/c. 

The master equation Q is valid even if the particle yields are initially not in full 
chemical equilibrium. It then dictates how the population changes with time. The BSW 
idea can now be formulated as follows: the overpopulation of meson states by a factor a, 
n n) K —> oi-n^ K , as compared to the standard equilibrium value in the hadron resonance 
gas model at T c should result in rates rescaled by a 5 , i.e., T^ od = a 5 ■ ^^k-^qb- ^ 
a = 2, the production rate for fls (as well as for other baryons and antibaryons) increases 
by a factor 32! One concludes, therefore, that baryons and antibaryons are readily 
produced by multi-meson reactions in an overpopulated bath of pions and kaons, and 
the equilibration should happen on timescales t eq ~ 1 fm/c |9:. The rapid (~ T 60 ) fall- 
off [Hj of the multi-particle reaction rates with temperature suggests that the "chemical 
freeze-out temperature" is also a measure of the phase transition temperature. 

However, there is a serious consequence of this argument: literally, it predicts 
that baryons and antibaryons are tremendously overpopulated as compared to their 
equilibrium values. Once produced, there is no way to get rid of them any more in 
the later stages, because the annihilation reactions are then suppressed dynamically. 
This is not seen in experiment. Inspect once again the master equation (JHJ): its (quasi-) 
stationary fixed point is reached when the expression in parentheses vanishes, i.e., when 
the annihilation rate is equal to the production rate. This corresponds to (anti-)baryon 
densities, rescaled by a factor (3, i.e., n s g — > (3 ■ n« B , where (3 = {af/ 2 . This may 
slightly change when other multi-meson channels are taken into account. For a = 2, 
(3 = 5.6, which predicts way too many (anti-)baryons in the system. In fact, the state 
with overpopulated (anti-)baryons is reached on a timescale of r = Tq/(3. This is short 
enough to compensate the rapid fireball expansion. However, once a large number of 
(anti-) baryons is produced, it is difficult to get rid of them quickly enough in order to 
reach standard hadron equilibrium values before the chemical freeze-out. The fraction 
of (anti-)baryons which annihilates is too small, because the corresponding annihilation 
rates are then not sufficiently large in an ordinary hadron resonance gas. 
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3. Fast equilibration due to Hagedorn States 

We propose to circumvent this drawbacks of the BSW scenario by postulating that 
the needed additional degrees of freedom close to T c are not light mesons ( "pions" and 
"kaons"), as proposed in Ref. [9j, but that they consist of heavy mesonic and baryonic 
resonances, called Hagedorn states. The conjectured unstable Hagedorn states can 
produce sufficiently many baryon-antibaryon pairs, on a sufficiently short time scale. HS 
efficiently and effectively can account for various multi-particle collisions ( "interactions" ) 
in a consistent way: 



This is akin to the Hagedorn's original idea (TO], that the strong interactions of low mass 
hadrons can be attributed to an exponentially increasing mass spectrum as T — > T c . 
This old idea has severe consequences for the time-scales of chemical equilibration - 
from the phase space arguments, one intuitively expects that multi-particle decays in 
(jUJ) dominate the BB production. 

Some comments to the Hagedorn picture are in order. Close to the crossover 
(critical) temperature, T c , the conventional description of QCD in terms of the known 
set of PDG hadron degrees of freedom can indeed not be sufficient (PDG: all known 
resonances from the particle data group booklet). In fact, one should insist that 
this set is not sufficient in order to understand the lattice results on the equation of 
state of QCD [HI El IT2] . The PDG-hadron resonance gas model falls short when 
estimating the pressure of lattice QCD as the crossover temperature is approached. 
This is not surprising, if one recalls the role of the highly lying massive states in the 
Hagedorn model. In the absence of a clear phase transition, it is easy to imagine 
that (additional, unknown) hadronic degrees of freedom can be used in the description 
of hot matter, in fact, even for temperatures above the crossover temperature |X3.. Of 
course, the further one goes away from "T c ", the less convenient such descriptions should 
become. Physically, this could be due to the increasing width of resonances [Ej. In the 
standard Hagedorn model, the density of states increases exponentially with their mass. 
Eventually, this leads to a divergence of the pressure, as the value of the temperature 
approaches the Hagedorn temperature from below. Is the Hagedorn model ruled out, 
therefore? We believe this is not necessarily the case. As pointed out in Ref. I14j. 
the model can be modified in a natural way that cures the problem with the pressure. In 
addition, quite recently the philosophy of Hagedorn behaviour has also become popular 
with regard to the equation of state of large N gauge theories [THj and Super- Yang-Mills 
theories [To] . 



HS «-> Til ■ 7T + Tl2 ■ K + 71% • K , 

HS <-+ B + B , 

HS «-> B + B + ni • 7T + n 2 ■ K + n 3 ■ K = B + B + X . 



(4) 
(5) 
(6) 
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3.1. Estimate of the baryon/antibaryon production 

The number density of HS states in the vicinity of the critical temperature is estimated as 
in Ref. [Oj: additional degrees of freedom to the PDG hadron resonance gas are needed in 
order to understand the energy density of a thermal system around T c , obtained in lattice 
QCD calculations. The additional energy density of order Ae c « 0.3 — 0.5 GeV/fm 3 is 
needed. Let us identify this with the energy density of the HS, Ae c = e HS . Than, a 
typical HS mass is M B .s ~ 3-6 GeV. The total number density of these HS is 

n HS «-£fi-« 0.05- 0.15 far 3 . (7) 

yyJ-HS) 

We now proceed to estimate the baryon-antibaryon production rate due to HS — > 
BB + X, where X stands for any possible number of additional hadrons. HS must 
be highly unstable: the phase space for multi-particle decays becomes immense with 
increasing mass. Extrapolating the width from the known meson resonances at 2 GeV 
(widths of 0.3-0.5 GeV) linearly (as suggested by the string model, e.g., in the last 
paper of Ref. [T^]), we find that the total width of such high mass mesonic states is 
~ 0-5 — 1 GeV. In the next subsection, we estimate the average baryon number (B) 
per unit decay of a HS within a microcanonical approach. Here we quote only the result: 
(B) « 0.2 — 0.4. Hence, the relative decay width is given by T BBX ~ (B) ■ m 100- 
300 MeV. Take the lower value in the following. Then, the production rate for baryon- 
antibaryon pairs is estimated 
dN, 



= T B°B d = ■ T HS ^ B§+X » 0.05 fin" 4 . (8) 

With n e ^ « n B « 0.04 far 3 at chemical freeze-out, at RHIC, for all baryons and 
anti-baryons, their chemical equilibration rate is given via r^g/nf « 1.25 fm" 1 : the 
chemical equilibration time is very small, T B h ? m « 0.8 fm/c. 

A large fraction of all possible reactions © also generates chemical equilibration 
among the "lighter" hadrons. For example, these distinct reaction channels easily alter 
the number of pions (or kaons) in the system (e.g., 7n <-> Air), so that in turn the light 
degrees of freedom become fully equilibrated chemically. 

The population of HS should die out very rapidly when the temperature is decreased 
within a narrow interval around the critical point. Thus, our scenario predicts that 
the chemical decoupling should happen naturally very close to T c , similarly as it was 
anticipated in the BSW scenario In contrast to the BSW mechanism, however, 
the approach proposed here does not lead to the BSW-oversaturation of baryons and 
antibaryons from decaying HS. Indeed, one might worry that such decays could produce 
too many BB pairs during the expansion and cooling of fireball. This is no problem 
as soon as the cooling proceeds smoothly, the densities of baryons and antibaryons stay 
close to their equilibrium values. However, in the extreme case of ultrafast cooling 
and decoupling, the number of additionally produced pairs (via the decay of all HS 
states) can be estimated to be Sn BB ~ (r ' bb+x/^hs) ' n HS ~ 0.2 nss- Thus, relative 
overpopulation is Sn BB /n e B ~ 0.2-0.4. As this is the most extreme case, the actual 
overpopulation will certainly be much smaller. 
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Figure 1. Multiplicity of individual hadrons coming per unit decay of a mesonic, 
nonstrange HS as a function of the mass E = Mhs- 



3.2. Microcanonical decay of a Hagedorn state 

We now provide an estimate for the individual branching ratios of the decay of a HS 
state into baryon-antibaryon pairs. The BB annihilation reactions at LEAR have been 
understood well within a statistical description [Tj\. Without any further information, 
we assume that the HS decay also with a statistical, microcanonical branching. 
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Figure 2. Branching probability distribution of additionally produced mesons HS — > 
B + B + n m ■ M as a function of the meson number per unit decay of a Mhs = 4 GeV 
Hagedorn state. 
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The microcanonical description of hadron reactions (e.g., pp or pp) developed in 
Ref. [T8| is used here. A state is described solely by its mass m = Mhs-, by its (reaction) 
volume V = Vhs = Mhs/c, where e denotes the mean energy density of a HS state, and 
by quantum numbers such as the net baryon number and the net strangeness. 

In Fig. ^ we depict the various numbers of hadrons h in the distinct decay channels 
HS — > h+X per unit decay. The HS states are purely mesonic, with no net strangeness. 
The energy density e of the HS states is varied in the range 0.25-0.75 GeV/fm 3 . The 
mass of the HS is varied continuously from Mhs — GeV to 10 GeV. For example, 
in a single decay of a state with mass runs — 4 GeV, there are about 2 tt + , 0.075- 
0.15 protons and neutrons each, 0.35 K + and 0.025-0.075 A's produced. For the mean 
baryon number, one finds (B) = 0.2-0.4. This was the input value used in the estimate 
of the partial decay width in the previous subsection. 

Fig.Elshows the branching distribution p{n m ) of the associated mesons as a function 
of the meson number per unit decay of a M^s — 4 GeV Hagedorn state. They are 
produced in HS — ► B + B + n m ■ M. The HS energy density is e = 0.5 GeV/fm 3 . The 
sum of the branching probabilities adds up to (B) « 0.3, as noted above. On average 
3 (stable) mesons accompany a HS decay into a baryon-antibaryon pair. Moreover, the 
direct decay, HS — * B + B, is not likely. 

3.3. Strange baryons and antibaryons 

The A particle has a sizeable production probability in the decay of nonstrange mesonic 
HS. One can estimate similarly that the A and the A can be quickly populated and reach 
their chemical equilibrium value by the interplay between the decay and the production 
of the HS. What about the most exotic baryonic state, the Q, with three units of 
strangeness? 

First, one might think that the Q can also be produced directly from nonstrange 
HS, i.e., via HS nonstrange ^ Q+B+X. For a M HS = 4 GeV HS, the average branching 
into fls per unit decay is calculated to be (Q) ~ 2 • 10 -4 , so TQ° d ^ 0.1 — 0.2 MeV. 
At T = 170 MeV, the density of fis in a PDG hadron gas is n£ « 4 • 10" 4 fm" 3 . The 
chemical equilibration rate is 

T chem = T prod ^ ^ 35 _ 5Q Me y ^ ^ ^ 4 _ g fm/c _ (g) 

n n 

This timescale is too long to (fully) explain a chemically saturated abundance. 

As a next step, consider (multi-) strange mesonic HS, for which the reaction 
HS(sssqqq) <-> Q + B + X might be sufficient. It turns out that (Q) w 0.05 for a 
m HS — 4 GeV and thus ^ 25 — 50 MeV. For the chemical equilibration rate one 
has 

T chem = TP ro d n HS (sssqqq) ^ 5Q _ ^ ^ Tn ^ 2 _ 4 fm/c (1Q) 
n n 

Here we have assumed that nHs(sssqqq)/nn nns nonstrange /Vb ~ 2.5. Altogether, 
strange Hagedorn states can explain the population of Q. 
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3.4- Importance of baryonic HS: from the Hagedorn temperature to the Hagedorn line? 

So far, we have only considered the role of mesonic HS. Now we discuss the role of 
baryonic HS. Until now, we have concentrated on the mid-rapidity region at RHIC, 
where //g ~ 0. Here, the fraction of baryonic to mesonic states can be estimated by 
comparison with that for normal hadrons: n bhs / n mhs ~ ^b/um ~ 0.04/0.3 = 0.13. 
One thus would conclude that these states do not carry a major fraction of the energy 
density. This can change, though, at finite \lb at AGS or SPS regime, where the baryonic 
HS achieve an extra enhancement factor e Ms / T . Such states could show up along a narrow 
//— T band below a critical "Hagedorn" line T H (fj, B ) [HI- This conjecture requires further 
development that is beyond the scope of this paper. 

Baryonic HS can also help to produce strange baryons. Consider the Q and any 
fi-like HS, respectively. The latter can decay and can be produced by the following two 
reactions: 

(a) QHS «-> B(^Q)+X (11) 

(b) QHS^Q + X (12) 

Here B ^ Q denotes baryon states with net strangeness less than three. The HS with 
higher net strangeness are populated by multi-hadronic fusion of baryons (resonances) 
with lesser strangeness and with other mesons, but with at least one kaon. On the 
other hand, the Qs then stem from the decay of these QHS. If such a scenario is to 
work, the QHS have to be produced sufficiently fastly. We now estimate the chemical 
equilibration times, both for saturating the QHS and for saturating the Qs. To do this 
one needs the relative branching probability separately for the decays happening via (a) 
and via (b). By employing the microcanonical model, one obtains for a Mqhs — 4 GeV 
Hagedorn state p( a ) = 0.9 and p^) = 0.1. The QHS predominantly decays into baryons 
with lesser strangeness. The chemical equilibration time for the QHS can be written 

(tuhZ ) ~* = FnHS-B+x « V{a) ■ T tot = 0.4 - 0.8 GeV . (13) 
Hence, such states do equilibrate very fast. For the Q we have 

I _chem\ 

[r n 1 = 1 qhs^q+x — -sr = P(b) • 1 — ■ i 14 ) 

The ratio of fi-like resonances to Q in the vicinity of T c is probably 0(1). We estimate 
ns/ijip + n n ) ps 4 at T=170 MeV. For nnHs/ n n = 1> one nas ~ 2 — 4 fm/c. Hence, 
exotic QHS can also be a key for explaining the chemical saturation of the Q. This 
would also be true, in particular, at SPS energies. Such a mechanism can work at finite 
[lb-, as the ratio nQHs/ n n does not depend on lib in Boltzmann approximation. 

4. Summary and Conclusions 

We have elaborated the special role of Hagedorn states for chemical equilibration of 
(strange) baryons and antibaryons. The chain of reactions © or © catalyzes rapid 
equilibration of antibaryons and baryons in the vicinity of the deconfinement transition. 
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The production and the decays of HS are governed by detailed balance, where both the 
continuous repopulation of HS as well as the annihilation of baryon-antibaryon pairs in 
the back reactions drives their chemical saturation. Three assumptions are necessary: 
(i) Ae HS « 0.3-0.5 GeV/fm 3 at T « T c ; (ii) r** s — 0.5-1 GeV; (iii) a microcanonical, 
statistical estimate of individual branching ratios. 

The Hagedorn states are additional degrees of freedom which represent complicated 
many-particle (hadronic or partonic) plasma correlations in hot dense matter. Close to 
the critical temperature this plasma is a strongly interacting phase, which contains such 
states. A clear cut proof, however, within the lattice QCD is not available at present 
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